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^ ■ Abstract 

We present the results of Suzaku X-ray follow-up observation of weak-lensing- 
detected halos in the field around galaxy cluster ZwC10823. 2+0425. We clearly 
detected X-ray emission associated with most of these halos and determined their 
detailed physical parameters such as X-ray luminosity, temperature, and metal abun- 
dance, for the first time. We find that the X-ray luminosity - temperature relation for 
these halos agrees with former typical results. With mass determined from the weak 
gravitational lensing data, the mass-temperature relation for them is also investigated 
and found to be consistent with the prediction from a simple self-similar model and 
results of the previous studies with both lensing and X-ray data. We would like to 
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emphasize that the self-similar scaling relation of mass and temperature is shown 
here for the first time using a weak-lensing selected sample, whereas previous stud- 
ies of the mass scaling relation used X-ray-selected samples of clusters. Therefore, 
our study demonstrates importance of X-ray follow-up observations of shear-selected 
clusters, and shows that a joint X-ray and lensing analysis will be crucial for clusters 
discovered by the forthcoming weak-lensing surveys, such as the one planned with 
Subaru/Hyper-Suprime-Cam. 

Key words: galaxies: clusters: individual (ZwC10823. 2+0425) — X-rays: galax- 
ies: clusters — gravitational lensing 

1. Introduction 

In the standard theory of structure formation in the cold dark matter (CDM) universe, 
larger structures such as dark halos corresponding to rich clusters of galaxies form through 
absorption of smaller dark halos and accretion of surrounding matter, which preferably occur 
along filaments associated with large scale structures in the universe. Although baryonic com- 
ponents such as diffuse hot gas and galaxies are visible in X-ray and optical bands, dark matter, 
which is a dominant component in mass, cannot be directly seen in any kind of radiation, and 
we recognized its presence only indirectly through dynamical properties of galaxies and hot gas. 

In recent years, however, this situation has been changed by the development of gravita- 
tional lensing techniques, which enable us to reveal matter distribution in the universe directly. 
We do not require any assumptions of the dynamical status and mass-luminosity ratio of the 
system to get the mass distribution (Schneider et al. 2006). Whereas strong lensing technique is 
essentially model dependent, weak gravitational lensing enables us to obtain mass distribution 
in the plane of the sky without any geometrical assumptions and measures directly galaxy clus- 
ter mass distribution out to the virial radius (Kaiser & Squires 1993; Bartelmann & Schneider 
2001). However, it is true that we cannot get detailed properties of baryonic components asso- 
ciated with dark matter halos by gravitational lensing alone. Therefore, joint studies of X-ray, 
optical and weak lensing datasets are very useful and provide us with important information 
of the interplay between baryonic and dark matter (Okabe & Umetsu 2008; Kawaharada et al. 
2010; Umetsu et al. 2010). 

In the very near future, large surveys in various wavebands (e.g. Subaru/Hyper-Suprime- 
Cam, eROSITA, SPT, and ACT) will provide us with huge datasets of galaxy clusters. In such 
future surveys, the weak gravitational lensing technique is one of the most powerful methods to 
find galaxy clusters, irrespective of the physical status of the baryonic components (Wittman et 
al. 2001; Miyazaki et al. 2002b; Miyazaki et al. 2007). Its detection efficiency of clusters/groups 
is a function of masses as well as the geometry of the universe. A typical procedure is as follows; 
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the lensing convergence field is reconstructed from galaxy's ellipticities convolved with a suit- 
able filter (Kaiser & Squires 1993; Seitz & Schneider 2001), and lensing signals above the noise 
of intrinsic ellipticities are identified as mass structures (Wittman et al. 2001; Miyazaki et al. 
2002b; Miyazaki et al. 2007). However, there are still some uncertainties in this method such 
as contaminations of member galaxies in shear catalogs, spurious peaks induced by intrinsic 
ellipticities, and the projection effects due to superpositions of structures at different redshifts 
(Miyazaki et al. 2002b; Hamana et al. 2009; Okabe et al. 2010a). Owing to not only these 
uncertainties characteristic of the weak lensing surveys but also the systematic difference of 
detection efficiency among weak lensing, X-ray, and Sunyaev-Zel'dovich effect (SZE) surveys, 
it is still unclear whether halo candidates detected by weak lensing always contain hot baryons. 
Thus, X-ray and/or SZE follow-up observations of halo candidates discovered by weak lens- 
ing analysis are very important as a pilot study before forthcoming large surveys in order to 
understand systematic differences among different wavelength surveys and the correlations of 
physical properties between hot baryons and total mass. In particular, a study of the scaling 
relations between lensing mass and baryonic observables is of prime importance for the cluster 
cosmology as well as understanding cluster baryonic evolution, such as radiative cooling, en- 
ergy feedback, and cluster mergers (Stanek et al. 2010; Vikhlinin et al. 2009a; Vikhlinin et al. 
2009b; Okabe et al. 2010c). 

We conducted Suzaku X-ray follow-up observation of weak-lensing-detected halo can- 
didates. X-ray Imaging Spectrometer (XIS) aboard Suzaku satellite (Mitsuda et al. 2007) is 
suitable for observing low surface brightness diffuse sources thanks to its low and stable back- 
ground (Koyama et al. 2007). Indeed, physical properties of the intracluster medium (ICM) 
around the virial radius, where X-ray emission is too faint for Chandra and XMM-Newton to 
get meaningful results, have been investigated for several clusters (Fujita et al. 2008; Bautz 
et al. 2009; George et al. 2009; Reiprich et al. 2009; Kawaharada et al. 2010; Hoshino et al. 
2010). This feature is also very useful to investigate diffuse X-ray emission associated with 
mass structures, even less massive objects, found by weak lensing surveys. 

As a pilot study for the above-mentioned purpose, we select a field around galaxy cluster 
ZwC10823. 2+0425. Several dark halos are found around the cluster by weak lensing study 
(Okabe et al. 2010b) as a collaboration of the Local Cluster Substructure Survey (LoCuSS; 
PI: Graham P. Smith). Figure 1 shows optical image of the field overlaid with the mass 
contours derived from the weak lensing analysis of Okabe et al. (2010b). Besides the central 
halo corresponding to the ZwC10823. 2+0425, four significant mass clumps are seen in the north, 
north-west, north-east, and south-east of it. We refer them as ZwC10823. 2+0425 (or C), N, 
NW, NE, and SE hereafter. The mass peak corresponding to Abell 664 is also found near the 
boundary. Two major red-sequence galaxy populations are found in this field, whose spatial 
distributions of both optical luminosity and number density are pretty similar to that of mass 
(Okabe et al. 2010b). The galaxies in the red-sequence at the lower redshift are apparently 
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associated with the C and NW halos. On the other hand, those in the higher redshift red- 
sequence seem to be related with the N and NE halos. The SDSS spectroscopic data are 
available for a few galaxies located in each halo center. We found that the mean redshift of 
galaxies associated with the C and NW halos is z = 0.2248, and that of N and NE halos is 
z = 0.472. The SE halo is likely to be associated with a bright galaxy at z = 0.103. However, a 
few galaxies in the background possibly contribute the lensing signal. This possible projection 
effect will be discussed in section 5. Although ZwC10823. 2+0425 is recognized as an X-ray 
source with the flux of 1.8 x 10~ 12 erg s" 1 cm" 2 in Bright Source Catalog of ROSAT All Sky 
Survey (Ebeling et al. 1998; Ebeling et al. 2000; Bohringer et al. 2004), its detailed spectral 
properties are not known. Moreover, X-ray emission associated with the other halos are not 
clearly detected. Therefore, deeper X-ray observation for this field is highly desirable to explore 
physical properties of the hot gas associated with the halos. 

In this paper, we present Suzaku X-ray follow-up observation of the field around 
ZwC10823. 2+0425 to investigate the physical properties of the hot gas in the dark matter 
halos found via gravitational lensing analysis. The rest of this paper is organized as follows. 
In section 2 we describe the observation and data reduction. In section 3 we present spectral 
analysis results. In section 4 we describe the imaging simulation to check the contamination 
of each halo's spectrum from the others. In section 5 we discuss the results and their implica- 
tions. In section 6 we summarize the results. Canonical cosmological parameters of H = 70 
Mpc -1 km s _1 , Q = 0.27, and A = 0.73 are used in this paper. Unless otherwise stated, all 
uncertainties are given at the 90% confidence level. 

2. Observation and Data Reduction 

We observed the field around ZwC10823.2+0425 with Suzaku on 2008 May 17-18. The 
field of view (FOV) of Suzaku XIS is shown in a Subaru optical image overlaid with the 
mass contours derived from the weak lensing analysis (Okabe et al. 2010b) in figure 1. The 
observation was performed at XIS nominal pointing. The XIS was operated in the normal 
full- frame clocking mode. The edit mode was 3x3 and 5x5, and we used combined data of 
both modes. The spaced-row charge injection was adopted for XIS. All data were processed 
with Suzaku pipeline processing, version 2.2. We employed calibration data files (20090925). 

The XIS data were processed through standard criteria as follows. Events with a GRADE 
of 0, 2, 3, 4, 6 and STATUS with 0:524287 were extracted. We excluded data obtained at the 
South Atlantic Anomaly (SAA), within 436s after the passage of SAA, and at low elevation 
angles from an Earth rim of < 5° and a Sun-lit Earth rim of < 20°. As a result, effective 
exposure time was 41.3 ks. Non X-ray background (NXB) spectra and images of XIS were 
generated using the ftool "xisnxbgen" (Tawa et al. 2008). Figure 2 represents an 0.5-8.0 keV 
XIS image combined from those of the front illuminated (FI) CCDs (XIS0, XIS3), overlaid 
with the mass contours. The image was corrected for exposure and vignetting effects after 
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Fig. 1. Optical image of the field around ZwC10823. 2+0425 overlaid with the mass contours derived from 
the weak lensing analysis. The field of view of Suzaku XIS is also shown in a black square. In addition to 
the central mass peak corresponding to ZwC10823. 2+0425, four mass peaks are seen around it. Hereafter, 
we refer them as C, N, NW, NE, and SE, as shown in this figure. The mass peak corresponding to Abell 
664 is also found outside of the XIS field of view. The contours are spaced in a unit of la reconstruction 
error, 5k = 0.02243. 

subtracting NXB, and smoothed by a Gaussian kernel with a = 0'.17. Enhanced X-ray emission 
from the ZwC10823. 2+0425, N and NE halo regions are clearly apparent in figure 2. 

3. Spectral Analysis 

For a spectral analysis of the XIS data, redistribution matrix files (RMFs) were generated 
using the ftool "xisrmfgen". In addition, ancillary response files (ARFs), which describe the 
response of X-Ray Telescope aboard Suzaku and the amount of the XIS optical blocking filters 
contamination, were generated with the ftool "xissimarfgen" (Ishisaki et al. 2007). XIS spectra 
for each sensor (XIS0, XIS1, XIS3) were fitted simultaneously For the spectral fitting, we used 
the energy band of 0.5 - 10.0 keV and 0.5 - 8.0 keV for FI CCDs (XIS0 and XIS3) and back 
illuminated (BI) one (XIS1), respectively. In the spectral fitting with Galactic absorption, we 
assumed Nh = 3.18 x 10 20 cm -2 (Dickey & Lockman 1990). The CXB level was estimated in the 
same way as Nakazawa et al. (2009) from the Lockman hole observation (Suzaku observation 
ID, 101002010). The detailed procedure is described in appendix 1 of Nakazawa et al. (2009). 
We defined the photon flux model as N(E) = 9.19 x 10~ 4 x E~ 1A in photons cm -2 s _1 keV -1 
FOV" 1 , where E is the photon energy in keV. 

3.1. The Background Model 

Before entering the spectral analysis of each halo, we need to construct a background 
model from astrophysical origin in addition to NXB. As shown in figure 3, we eliminated regions 
corresponding to each halo and the foreground bright star near the west edge of the XIS FOV, 
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Fig. 2. XIS image of the field around ZwC10823.2+0425 in the 0.5-8.0 kcV band combined from XIS 
FI CCD images, overlaid with the mass contours. The image was corrected for exposure and vignetting 
effects after subtracting NXB, and smoothed by a Gaussian kernel with a = 0'17. 

the rest of which was used for spectral fit to determine the background model. For the analysis 
of the background model, uniform emission over a circular region with 20' radius was used as 
an input image to generate an ARF. The Cosmic X-ray background (CXB), the Galactic halo's 
hot gas (GH), and the local hot bubble (LHB) were considered as the background components. 
Then, the spectrum of the background region was fitted by a model as follows, 

aped + wabs x [apecl + apec3 + power law) (1) 

where apecl, apecl + apec3, and powerlaw represent the LHB, GH, and CXB, respectively. 
The temperature of the LHB was fixed to be 0.08 keV. The metal abundance of the both LHB 
and GH were also fixed to be solar. 

The spectrum of the background region fitted with the above-mentioned model is shown 
in figure 4, where the black, red, and green crosses show the spectrum of XIS0, XIS1, and XIS3, 
respectively. The total and each component of the best fit model spectra are also plotted as 
solid and dashed histograms, respectively. The detailed results of the fit are summarized in 
table 1. At first glance, the higher temperature component of GH may seem to be too high. 
We tried to fit the background region spectrum by a model with a single temperature GH. 
However, the results were not so good and we had some residuals around ~ 2 keV in the data. 
This might be because our spectral modeling is too simple for the data. However, our main 
purpose here is to construct a plausible background spectrum model to investigate the halos' 
hot gas. Thus, we do not pursue this issue in more detail. 
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Fig. 3. A Region used to determine the background model (yellow field). We eliminated regions corre- 
sponding to each halo and the foreground bright star near the west edge of the XIS FOV (green circles 
and ellipses). 





m 




iJiilJ 1 iiij 

ij J i J 




iff 





1 2 5 

channel energy (keV) 



Fig. 4. XIS spectrum of the region shown in hgurc 3 fitted with the background model described in the 
text. The black, red, and green crosses show the spectrum of XISO, XIS1, and XIS3, respectively. The 
total and each component of the best fit model spectra are also plotted as solid and dashed histograms, 
respectively. 
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Table 1. Best fit parameters for the XIS spectrum of the background region. 



The background model 





0.08(fixcd) 






i.i9±g:|g x 10 


-2 




S4+ - 24 




A^GH.low 


5.76^:24 X 10 


-4 


fc^GH.high* 


i 71 +0.45 
L - ' -0.38 




-^GH.high* 


5.18±^j? x 10 


-4 


Tpl* 


1.4(fixed) 




A^PL § 


9.19 x 10- 4 (fixed) 


X 2 /d.oi. 


403.5/336 





* Temperature of the each component in keV. 

t Normalization in the apec code for each component. 

' Photon index of the power-law component. 

§ Normalization in the power-law component. 

3.2. Hot Gas in Each Halo 

We performed spectral analysis for each halo using the background model in the previous 
subsection. Figure 5 shows regions used in the analysis of each halos, where the center of each 
region was determined through the each halo's mass peak in the weak lensing data of Okabe 
et al. (2010b). The radii of the regions are 2.5', 2.0', 1.8', 1.8', and 2.0' for C, N, NE, NW, 
and SE, respectively. It is true that these radii should be as large as the expected virial radii 
of the corresponding halo. However, because each halo partially overlaps with one another 
and the spatial resolution of Suzaku is moderate, the radii have to be small enough to reduce 
contamination from the others. No core removal was made in the spectral extraction. We 
adopted model images of the /3-model profile with core radius 0.2 Mpc and /3 = 0.66 at each 
halo's redshift as an input image to generate ARFs, considering that Suzaku does not have 
spatial resolution good enough to resolve inner structures of each halo and that there is no 
X-ray image of this field deep enough for our purpose by other instruments. 

Systematic errors of both CXB and NXB were taken into account in the following anal- 
ysis. It is well-known that the CXB fluctuations can be modeled as ctcxb/^cxb oc fl~ 0S S^ 25 , 
where Q e and S c are the effective solid angle and upper cutoff flux of a point source, respec- 
tively From the HEAO-1 A2 results, ct C xb//cxb = 2.8% with Q e = 15.8 deg 2 and S c = 8 x 10" 11 
erg s -1 cm -2 (Shafer 1983). We adopted upper cutoff flux S c ~ 1.0 x 10~ 14 erg s -1 cm -2 taking 
into account the fact that point sources with similar flux were clearly detected in the Lockman 
hole observation (Suzaku observation ID, 101002010). Thus, the CXB fluctuations for C, N, 
NE, NW are expected to be 26, 33, 36, and 36 % at the 90% confidence level. Tawa et al. 
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Fig. 5. Regions used in the spectral analysis of each halo. 

(2008) reported that reproducibility of NXB was 6.0% and 12.5% for XIS FI and BI at the 90% 
confidence level, respectively. We did not take account of systematic errors of LHB and GH, 
considering that Yoshino et al. (2009) reported that prominent spatial fluctuations in a scale 
smaller than that of XIS FOV were hardly seen in LHB and GH. In the spectrum analysis, 
NXB and CXB components were fluctuated at the 90 % confidence level of the systematic 
uncertainty mentioned above, which caused the changes of the best-fit parameters in the fits 
and gave us the systematic errors. 

Each halo's spectrum was fitted by the photoabsorbed single temperature apec model 
(wabs x apec) with the background model described in the previous subsection. The redshift of 
each halo was fixed to be a value obtained from the SDSS data. We fixed the metal abundance 
for NE, NW, and SW to be 0.3 solar because we do not have enough photons for these regions 
to determine it. Figure 6, 7, 8, 9, and 10 represent the XIS spectrum fitted with the above- 
mentioned model of C, N, NE, NW, and SE, respectively. In each figure, the black, red, and 
green crosses show the spectrum of XIS0, XIS1, and XIS3, respectively. The best fit model 
spectra are also plotted as solid histograms. In general, all data except for SE are fitted well 
by the adopted model. Unfortunately, however, the statics of the SE spectrum is not good 
enough to obtain meaningful fitting results. This is consistent with the fact that SE lensing 
signal in the mass map is due to the projection effect, which is mainly associated with a galaxy 
at low redshift z ~ 0.1, as described in section 5. The detailed fitting results are presented in 
table 2. The luminosity in 0.5 — 10.0 keV of each halo is calculated from the fitting results and 
summarized in table 3. 

We used model images to generate ARFs because Suzaku has only moderate spatial 
resolution (~ 2'), and because we do not have appropriate X-ray images with better spatial 
resolution by other instrument for this field. However, it should be noted that an effective area 
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Fig. 6. XIS spectrum of the ZwC10823. 2+0425 fitted with the model described in the text. The black, 
red, and green crosses show the spectrum of XISO, XIS1, and XIS3, respectively. The best fit model 
spectra are plotted as solid histograms. 




Fig. 7. Same as figure 6, but for the N halo. 
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Fig. 8. Same as figure 6, but for the NE halo. 
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Fig. 9. Same as figure 6, but for the NW halo. 
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Fig. 10. Same as figure 6, but for the SE halo. Unfortunately, however, the statics of the SE spectrum is 
not good enough to obtain meaningful fitting results. 



Table 2. Best fit parameters for the XIS spectrum of each halo region presented in figure 6, 7, 8, and 9. 
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o 7Q+1.26+0.S6 
owy_ .go_o.66 

0.3(fixed) 
0.472(fixcd) 
4.04 



0.50+0.22 1f) - 
0.56-0.35 A 1U 



48.9/60.0 



n , +0.75+0.57 
z - u ^-0.44-0.42 

0.3(fixed) 
0.2248(fixcd) 

1 7 o+0.23+0.37 1f) - 
i - ' J -0.36-0.46 A iu 



40.8/52 



* Temperature of the each halo in keV. Errors are 90% statistical and systematic, respectively. The latter is evaluated by shifting the NXB 
and CXB within their systematic errors, 
t Metal abundance. 

t Redshift obtained from the SDSS data for each halo. 
§ Normalization in the apec code. 
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Table 3. Luminosity in 0.5-10.0 keV of each halo derived from the fitting results listed in table 2 with statistical and 
systematic errors. 



C N NE NW 

l,™;T,ndhJm 44 Bro s" 1 "! A 7fi+°- 11+0-08 OOO+0.44+0.34 1 OQ+0.21+0.25 n 1 c+0.05+0.' 

Luminosity (IU ergs j 4.7b_ 14 _ 08 »-«^_ .57-o.3i Lm -o. 25-0. 33 u - i5 -o.05-o. 



for diffuse sources depends on their spatial distribution. Therefore, absolute values of the apec 
code's normalizations in table 2 and fluxes listed in table 3 could change if the adopted model 
images change. To evaluate this, we generated ARFs with a point source and /3-model images 
with different core radii (r c = 0.1 and 0.3 Mpc). We confirmed that this changed the fluxes by 
~ 50 % at most. Thus, we believe that uncertainty due to image modeling for ARF generation 
is at a few tens of percent level. 

Some of the halos might have central cool cores that are not resolved in our observation. 
This means that the halos' temperature and metal abundance are possibly underestimated and 
overestimated, respectively. Furthermore, unresolved point sources with harder spectrum such 
as AGNs could be contaminated, which results in overestimation of both the ICM temperature 
and metal abundance. If these kinds of unresolved contamination sources are dominant, the 
observed metal abundance could be higher than that of a typical cluster without a cooling core. 
However, such a trend is not seen at least in the central cluster and N halo. Thus, we safely 
conclude that the effects are quite limited at least for these halos. 

4. Imaging Simulation 

Because each halo might partially overlap with one another, we used regions with radii 
of only ~ 2' as shown in figure 5. Still, each halo's spectrum could be contaminated by the 
other halos and point source near the west edge of XIS FOV. We simulated X-ray emission 
from them with the ftool "xissim" in order to estimate contamination to each halo region from 
the others. We placed five /3-model X-ray halos with r c = 0.2 Mpc at the position of C, N, NE, 
NW, and SE, and one point source at the position of the bright star. As a spectrum model, we 
used the best-fit models listed in table 2 for C, N, NE, and NW. We also fit the SE and point 
source in a similar way, and the best-fit results were used as a spectrum model for the image 
simulation. The resultant simulated image of XIS0 is shown in figure 11. 

We checked contamination to each halo region from the other halos and point sources. 
We confirmed that only 4 % of photons detected in the C region was from the others. For N, 
16 % of photons were from the others, most of which was originated from C. 32 % of photons 
detected in NE regions were from the others, about two thirds and one third of which were from 
C and N, respectively. As for NW, 56 % of photons detected in the region was from others, 
about three quarters and one quarter of which were originated from C and N, respectively. 
Even in case of the faintest NW halo, contamination from the others alone cannot account 
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Fig. 11. Simulated image of XISO. 



for the detected X-ray emission. It should be noted that relatively faint halos' (NE and NW) 
spectrum are significantly contaminated by those of the brighter ones (C and N) with higher 
temperature, which means that the temperatures of NE and NW are probably overestimated. 
This effect is not corrected in the spectrum analysis because our imaging simulation is just aimed 
at confirming that the fainter halo's emission cannot be accounted for only the contamination 
from the brighter ones and that the hot gas associated with these fainter halos is surely detected. 

5. Discussion 

First, we investigate the X-ray luminosity-temperature {Lx — k-QT) relation for the halos. 
Figure 12 shows the Lx — k^T relation obtained from our results, where the errors are converted 
into at the one-sigma level. The best-fit power-law function is 

logL x (10 44 erg s" 1 )^)- 1 = -1.59+°f + 3.21+^fog fc B T(keV), (2) 

where 

E(z) = y/n (l + z)* + A . (3) 

Although errors are large because of relatively poor statistics, our results are consistent with 
the former typical ones, which show LxE(z)~ l oc T 2 ' 5-3 (e.g., Ikebe et al. 2002; Branchesi et 
al. 2007; O'Hara et al. 2007). 

Second, let us investigate the mass-temperature (M — k^T) relation. Here, we made a 
secure selection of background shear catalog in a color magnitude plane in order to avoid con- 
taminations of member galaxies associated with each clump (Umetsu & Broadhurst 2008; Okabe 
et al. 2010b). As described in section 1, there are two red-sequences at z ~ 0.22 and ~ 0.47. 
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Fig. 12. X-ray luminosity - temperature relation for the halos. Errors are converted into the one-sigma 
level. The best- fit power-law function (equation (2)) is also overlaid. 

Thus, we conservatively cut the broad range of the color. The resultant number density of the 
background galaxies decreased to one- fourth of the original shear catalog [n g ~ 7.6 [arcmin -2 ]). 
We calibrated the mean redshifts based on the COSMOS photometric redshift catalog (Ilbert 
et al. 2009). The existence of halos at various redshifts makes it difficult to conduct one- 
dimensional tangential distortion analysis in order to measure each halo mass. In contrast, 
using two-dimensional lensing distortion pattern, we are able to model lensing signals from 
all halos as well as the projection effect. Therefore, we conducted multi-components fitting of 
the two-dimensional lensing distortion pattern. Here, we briefly summarize the fitting method. 
Full descriptions will be presented in the forthcoming paper (Okabe et al., in preparation). 
We first pixelized the shear distortion pattern into a regular grid of 2'. x 2'. without any spatial 
smoothing procedure, considering the limited number of the available background galaxies. The 
pixelized distortion signals, (g a ), and the mean position, (9 a ), in the n-th grid are calculated 
with statistical weight, Ui, 



v) 



E UiQoc,i 

E ^ 
E ^ 

1 



U 



a 



9,i 



(4) 



(5) 



(6) 



where g a>i and 9 a i are the reduced shear and position of the i-th galaxy in the n-th grid, 
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respectively. a g ^ is the rms error for shear measurement, and a g = {^cr^^/N) 1 ' 2 ~ 0.48 is 
the softening constant variance (Okabe & Umetsu 2008; Okabe et al. 2010b). Note that this 
pixelization method is different from that of making mass maps such as figure 1, where we 
pixelized the shear distortion pattern with a Gaussian smoothing kernel (FWHM= 1^67) using 
a shear catalog without the color selection (Okabe et al. 2010b). We then conducted the \ 2 
fitting given by 

2 Npixel 

X 2 = E E (^(^n) -^(0n;p) (m ° del) )^(^)(^(^) -^(0n;p) (mOdeI) ), (7) 

where p is the parameter vector. C a p is the error covariant matrix of shape measurements in 
the form of C a p(O n ) = S^pag(O n ), where S^p and (J 2 (0 n ) are a Kronecker delta function and the 
statistical error of the pixelized shear (Oguri et al. 2010), respectively. 

We considered all mass components of C, N, NE, NW, SE and A664. We adopted an 
NFW mass profile (Navarro, Frenk, & White 1996; Navarro, Frenk, & White 1997) for each 
halo model. The NFW profile can be described with the two parameters; the virial mass M vir 
and concentration parameter c v - 1T . According to cosmological A r -body simulations based on 
the CDM universe, it is well-known that the concentration parameter is a weakly decreasing 
function of the virial mass M vir and redshift z, e.g., c vir = 7.85(M vir /2 x 1O 12 M )~ ' 081 (l + z)~ - 71 
(Duffy et al. 2008). We here assumed this relation for halos of N, NE, NW, SE and A664. This 
is because their angular sizes are likely too small to constrain the concentration parameters well, 
taking into account the fact that they are at higher redshift (z ~ 0.4) or less massive objects. 
On the other hand, we treated the concentration parameter as a free parameter for the nearby 
main cluster, C. In addition, the positions of the halo's center (x c ,y c ) for C, N, NE, and NW 
were also parameterized so as to avoid a mis-center in mass measurements. The center positions 
for the other components were fixed at those of the brightest galaxy appeared in the optical 
image. As a result, we used 15 parameters to model the two-dimensional weak-lensing analysis 
data. The modeling shear distortion pattern, g a (^n]p)^ modcl \ is described by a summation of 
shear signals of each mass component. We adopted the Markov Chain Monte Carlo (MCMC) 
method with standard Metropolis-Hastings sampling. We restricted the sampling range to 
M vir < 5 x lO 15 /i _1 M , c vir < 50, |x c -x peak | < V. and |y c -2/ p0 ak| < 1', where x pcak and y pcak are 
the peak coordinates appeared in weak-lensing mass map. 

The obtained virial masses M vir and M 50 o, mass within a sphere of the radius r 50 o are 
shown in table 4, where r^o is the radius within which the mean density is 500 times of the 
critical density of the universe. The parameters of NW and SE halos were not constrained 
well by this analysis, because their lensing signals are much smaller than those of the other 
massive objects. The mode of the posterior distribution for the SE halo's virial mass is M V - 1T ~ 
6 x 10 13 /i -1 M Q . This is typical to galaxy groups rather than clusters, which implies that it 
would be mainly due to the lensing signal from the mass structure associated with the galaxy 
at z — 0.103 and possible background galaxies. To investigate this possible projection effect 
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Table 4. The best-fit masses (M V [ T and M500) for clump candidates in the unit of 10 14 h 1 Mq. 



Halo Name 


C 


N 


NE 


NW 


SE 


A664 


M vir 
M500 


n QO+0.90 

o-oo_ 1 47 
2 40+ - 64 


c- 0O+2.24 

°- oz -3.07 

q 91 +1.19 

°' zi -1.66 


9 01 +1.19 
z - 01 -1.98 

1 01 +0.65 
1 ' 0i -l.ll 


< 0.43 

< 0.25 


< 1.05 

< 0.58 


Q 1 n + 2.72 

8 - iz -3.29 

A Q9 + 1 - 36 

*' JZ -1.75 



in the SE clump region, we calculated number densities of galaxies in a given photometric 
redshift slice with SDSS photometric data in the same way of Kawaharada et al. (2010), taking 
account of an uncertainty of the photometric redshifts. We calculated peak significance values 
based on the mean density in the Subaru's field-of-view, in order to avoid redshift evolution 
of the number density of galaxies and physical scale difference. We found a ~ 3.5cr peak at 
z = 0.1 with a slice of 8z = (500 kms -1 )/ (1 + z)/c ~ 0.002. We also found another ~ 3a peak 
in the redshift range of \z — 0.4 1 < 5z = 0.01, which corresponds to the velocity dispersion of 
~ 2000 kms -1 . Therefore, SE clump is likely to be a superposition of mass structures at these 
two different redshifts. Since a part of the NW halo overlaps with the main cluster C within 
rsoo, the NW halo might be a tidally-truncated substructure. However, we cannot rule out 
a possibility that all mass halos are isolated objects because we cannot distinguish structures 
along the line-of-sight in a few tens of Mpc scale observationally. 

Figure 13 shows the M — ksT relation for the halos. Again, the errors are at the one- 
sigma level. It should be noted that the mass used here is obtained from the weak-lensing data 
alone, which is perfectly independent of the X-ray results and free from the assumption about 
the dynamical status of the system. First, we fit the data into a power-law function setting 
both the index and normalization as a free parameter. The best fit power-law function is 



2.38 



-0.78 
-0.95 



M 500 E(z) = 2.63l°J* (J^J x lO^Mo, (8) 

which is plotted as solid lines in figure 13. Here, we do not take into account intrinsic scatter 
of the scaling relation (Okabe et al. 2010c). This result is consistent within large errors with 
what a self-similar model predicts, ME(z) oc T 15 (Kaiser 1986). Next, we fit the data into a 
power-law function with the index fixed to that of the self-similar model. The resultant function 
is 

M 500 E(z) = 2.47i°$ (^) 3/2 x lO^Mo, (9) 

which is also plotted as dashed lines in figure 13. Interestingly, the obtained normalization 
is consistent with that of the previous study (2.45^o;24;Okabe et al. 2010c) for a sample of 12 
clusters based on weak-lensing and X-ray datasets, which is plotted as dotted lines in figure 13, 
though they almost overlap with the dashed lines. Although we assume the mass-concentration 
relation for three halos, note that the masses at the intermediate scale such as r 50 o hardly depend 
on the assumed relation. As mentioned in section 3.2, the possible presence of cool cores and/or 
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k B r [keV] 

Fig. 13. Mass - temperature relation for the halos. The errors are at the one-sigma level. The best-fit 
power-law functions (equations (8) and (9)) are also overlaid with solid and dashed lines, respectively. 
The dotted lines denote the result of Okabe et al. (2010c) for a sample of 12 clusters based on a joint 
weak-lensing and X-ray analysis, which almost overlap with the dashed lines. 

AGNs has very limited impacts on the temperature measurement. We would like to emphasize 
that this is the first confirmation of a self-similar scaling relation of mass and temperature 
for halos discovered by weak-lensing signals, whereas samples of X-ray selected clusters were 
used in previous studies of the mass scaling relation. The NW halo's mass expected from its 
temperature via equation (9) is twice as large as its actual upper limit, which might indicate that 
the NW halo is tidally truncated by the primary halo, ZwC10823. 2+0425. Another possibility 
is the probable overestimation of the temperature because of the contamination as mentioned 
in section 4. 

Halo detection criteria in weak lensing survey might be affected by the Gaussian smooth- 
ing scale adopted in mass reconstruction. We therefore investigated dependence of significance 
level of peaks on smoothing scales. Here, we do not make a selection of background galaxies, 
as a usual way. The significance levels changed by only < 15% as the smoothing scale changed 
by ~ 30%, and thus we conclude that this effect is negligible for halo detections. 

In order to compared the significant level of peaks between the results with and without 
the color selection, we made the projected mass density map for the best-fit mass model with 
the same smoothing kernel as figure 1. This is because the number density of the secure 
background galaxies is too small to reconstruct the mass map with the smoothing scale same 
as figure 1. The peak k for C and NE halos in the model map are higher than by 2.6a and 1.5<r 
level of the reconstruction errors (5k ~ 0.022) of figure 1. This implies that a contamination of 
the member galaxies in the shear catalog without color selection (known as a dilution effect) is 
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significant for C and NE halos. The NW peak in the model map is about la level lower, while 
SE and N halo peaks are comparable between two maps. Although we do not have a definite 
idea why the significance of NW halo decreases, it is possible that the background selection in 
the color-magnitude plane did not work well for some reasons. To make more secular selection 
of background galaxies, we need more detailed multi-band imaging data, as demonstrated by 
Umetsu et al. (2010). We also computed the contribution for the lensing signal from each halo 
to investigate the projection effect. The lensing effect from the other halos on C and N halo 
regions is only about 10%. In contrast, the contribution from the others accounts for ~ 20%, 
~ 20% and ~ 40% for SE, NE and NW halo regions, respectively. As for the observed lensing 
signals of these three halo regions, therefore, the contribution from the others is not negligible. 
In particular, the projection effect and/or contribution from nearby halos are critical for halo 
detection with the weak-lensing technique, because they boost the peak strength in the mass 
map. For instance, the peak projected density of NW halo alone is only at ~ 2a level of the 
reconstruction error of figure 1. This indicates that we need further careful lensing studies for 
halo candidates discovered by lensing maps. Otherwise, we cannot rule out possibility that 
shear-selected catalog of clusters/groups is contaminated with some less massive objects. As 
for cluster surveys with weak lensing technique alone, this difficulty should be kept in mind to 
achieve a complete sample. 

6. Conclusions 

We observed the field around the galaxy cluster ZwC10823. 2+0425, where several dark 
matter halos were found with the weak gravitational lensing technique. X-ray emission asso- 
ciated with most of these halos was clearly detected. From spectral analysis of the XIS data, 
we determined temperature for the all halos detected in X-ray, and metal abundance for some 
of them, for the first time. With these results, we obtained the X-ray luminosity-temperature 
relation, which is consistent with the former typical results. The gravitational mass of each 
halo was obtained from the weak lensing data, which is absolutely independent from the X-ray 
data and free from the assumptions about the dynamical status of the systems. With the 
lensing and X-ray results, we investigated the mass-temperature relation for the halos, which is 
found to be consistent with the prediction of a simple self-similar model and previous studies. 
We first found the self-similar scaling relation of temperature and mass for halos detected by 
weak-lensing signals. Therefore, our study demonstrates the importance of X-ray follow-up 
observations of shear-selected clusters, and shows that a joint X-ray and lensing analysis is cru- 
cial for clusters discovered by the forthcoming weak-lensing surveys, such as the one planned 
with Subaru/Hyper-Suprime-Cam. It is expected that the weak-lensing surveys will provide 
us with the uniform catalog of galaxy clusters irrespective of cluster dynamical states, though 
there are still some difficulties as discussed in section 5. This is because the lensing technique 
depends on the cluster masses alone, which is free from the dynamical state. In contrast, X-ray 
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and SZE observables are obviously affected by cluster merger phenomena (Ricker & Sarazin 
2001; Takizawa et al. 2010). A study of mass-scaling relations with statistically complete sam- 
ples of galaxy clusters will enable us to understand cluster baryonic evolution, such as radiative 
cooling, pre-heating, and cluster mergers (Stanek et al. 2010), and to make a robust mass-proxy 
based on principal component analysis (Okabe et al. 2010c). This will be an important step 
to constrain the equation of state for the dark energy through the cluster mass function. The 
future X-ray projects, such as Astro-H and eROSITA, will be powerful in this regard. 
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